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ABSTRACT 

We present time-series data on Jupiter Family Comets (JFCs) 17P/Holmcs, 
47P/ Ashbrook- Jackson and 1 3 7P/ Shoemaker-Levy 2. In addition we also present re- 
sults from 'snap-shot' observations of comets 43P/ Wolf-Harrington, 44P/Reinmuth 2, 
103P/Hartley 2 and 104P/Kowal 2 taken during the same run. The comets were at 
heliocentric distances of between 3 and 7 AU at this time. We present measurements 
of size and activity levels for the snap-shot targets. The time-series data allow us to 
constrain rotation periods and shapes, and thus bulk densities. We also measure colour 
indices (V — R) and (R — I) and reliable radii for these comets. We compare all of 
our findings to date with similar results for other comets and Kuiper Belt Objects 
(KBOs). We find that the rotational properties of nuclei and KBOs are very similar, 
that there is evidence for a cut-off in bulk densities at ~ 0.6 g cm -3 in both popu- 
lations, and the colours of the two populations show similar correlations. For JFCs 
there is no observational evidence for the optical colours being dependant on either 
position in the orbit or on orbital parameters. 

Key words: comets: general - comets: individual: 17P/Holmes, 47P/Ashbrook- 
Jackson, 1 3 7P/ Shoemaker-Levy 2 - Kuiper Belt - techniques: photometric 



1 INTRODUCTION 

It was only in the 1950s that the modern description of 
comets, as bodies with a small nucleus of ice and rock that 
produce a coma through sublimation, and tails by solar win d 
pressure on this coma, was postulated JWhippleil95dll95lD . 
Furthermore, it is difficult to study the nucleus itself, as 
comets are generally active (outgassing to produce a coma) 
when near to the Earth, and when active the flux reflected 
by the coma greatly exceeds the flux from the nucleus. 

There are three techniques by which this difficulty 
can be overcome at optical wavelengths. First, the nu- 
cleus can be seen within the coma by employing space- 
craft to image the comet from very close range. This was 
achieved with great success by the VEG A and Giotto 
missi o ns at comet lP/H alley during 1986 dSaedeev et alJ 
Il98rj; iKeller et alJ Il986l) . proving Whipple's 'dirty snow- 
ball' model, and by the recent mis sions Deep Space 1 at 
19 P/Borrellv JSoderblom et all20o5) . Stardust at 81P/Wild 

2 ferownlee et al.ll2004l) and Deev Imvact at 9P/Temple 
1, the last of these making contact with the nucleus 
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in spectacular fashion jA'Hearn et al.lf2005l) . The Rosetta 
space-craft is on its way to 67P/Churyumov-Gerasimenko 
iSchwehm fc Sc hulz 1998^), and will map its nucleus in de- 
tail and land a probe on the surface. While these missions 
provide the most complete information on the nucleus, they 
are expensive and limited in the number of targets they can 
observe. 

A second method relies not on going to the comet, but 
waiting for the comet to come to us. By taking advantage of 
very close passes of comets to the Earth, and very high reso- 
lution imaging using the Hubble Space Telescope (HST), the 
signal from th e nucleus can be se parated from the surround- 
ing coma (e.g. lLamv et alJll999l) . This method can only be 
applied to a limited number of comets, due to the rarity of 
close passes of comets to the Earth. 

The third method, which we employ here, and have pre- 
viously applie d to comets 7P/Pons-Winnecke, 14P/Wo lf and 
92P/Sanguin iSnodgrass. Fitzsimmons fc Lowrvl2 005'). is to 
image comets when inactive at large distances. Our pro- 
gramme of obtaining nuclear properties has concentrated 
on Jupiter Family Comets (JFCs), which are short pe- 
riod comets whose orbital evolution is dominated by the 
planet Jupiter. Believ ed to come from the Kuiper Belt 
illp fc FernandezHl997T) . these nuclei are left over remnants 
from the formation of the solar system. JFCs are generally 
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found to be active when at heliocentric distances of < 3 
AU, as within this distance the sublimation rate for water 
ice (the primary volatile constituent of nuclei) increases dra- 
matically. This method is difficult due to the small size (radii 
of a few kilometres) and very low albedos (typically ~ 4%) 
of nuclei, which makes them very faint targets at large dis- 
tances. However, by careful selection of targets with good 
pr evious absolute ma gnitude estimates (i.e. the compilation 
bv lLamv et al.lf2004l) that have > 3 AU, we can find 
inactive nuclei which are observable at sufficient S/N using 
4m class telescopes. 

We employ time-series photometry to constrain the 
shape and bulk density of nuclei; taking data over a series of 
days to find the rotation period from the periodic variations 
in the brightness of the nucleus. Additionally, 'snap-shot' 
observations only give an effective radius at one generally 
unknown rotational phase. The time-series data allows the 
true effective radius to be measured. Here we present time- 
series photometry of three JFC nuclei, in addition to results 
from 'snap-shot' observations of four others. 

In section [5] we describe the observations and data re- 
duction procedures employed in this work. We explain our 
analysis methods, and present our results from the snap-shot 
and time-series data, in sections [3] and 2] Section|^|discusses 
our findings in terms of results for other comets, draws con- 
clusions about the population of JFCs as a whole, and com- 
pares JFCs with Kuiper Belt Objects (KBOs), the supposed 
parent population for JFCs. Finally, sectionE]summarises all 
of our results. 



2 OBSERVATIONS AND DATA REDUCTION 

We used the 3.6m New Technology Telescope (NTT) at the 
European Southern Observatory's La Silla site on the nights 
of the 5th to 7th March 2005. The later half of the 1st night 
and the third night were lost to poor weather, but conditions 
were photometric during the remaining time. The comets 
targeted were at heliocentric distances between 3 and 7 AU; 
tableQgives the orbital positions of the comets at this time 
and a summary of the exposures taken. CCD imaging was 
performed using the red arm of the EMMI instrument which 
is mounted at the f/11 Naysmith-B focus of the NTT. The 
EMMI red arm contains a mosaic of 2 MIT/LL 2048x4096 
CCDs, and was used in 2x2 binning mode to give a pixel 
scale of 0.332 arcsec per pixel. The effective field of view is 
9.1x9.9 arcmin 2 . The images were taken primarily through 
the Bessell R band filter, with at least one set of images 
taken through the Bessell V and / filters each night for each 
comet to allow measurement of colour indices. All images 
were taken with the telescope tracking at the sidereal rate, 
with exposure times chosen so that the apparent motion of 
the comet would be less than 0.5", and would thus remain 
within the seeing disk. This allowed us to perform accurate 
measurement of the stellar-background point spread func- 
tion (PSF), and highly accurate differential photometry. The 
FWHM of the stellar-background PSF was found to vary 
between 0.6" and 3.1" over the two nights, with a median 
of 1.1". The U_R/-filter exposure times used are listed in ta- 
ble Q We observed each comet in blocks of exposures lasting 
~ 20 minutes, and cycled between each of the comets that 
were visible at any given point in the night. This procedure 



gave us good temporal coverage without losing too much 
time slewing the telescope between targets. 

T he reduction was performed using standard IRAF 
tasks IITodvl Il986l . Il993l) . Bias-subtraction was performed 
utilising the prescan strip of CCD 1. Twilight sky exposures 
were acquired during the evening and morning, and these 
were combined to produce a flat-field imag e for each filter 
for each night. We have previously JSnodgrass et alJl2005h 
found that there is a variation in illumination of NTT flat- 
field images with rotator angle, due to scattered light from 
the Naysmith mirror baffle. To allow for this, we took each 
set of flats at a number of rotator angles, which when aver- 
aged correct the problem. 

Differential photometry was performed on the comets 
using the IRAF packages DIGIPHOT and APPHOT dPavisl 
1999). Aperture radii were set equal (to the nearest integer 
pixel) to the FWHM of the stellar PSF, which has previ- 
ously been fo und to be the optimum for maximising S/N 
llHowel]|ll989D . Choosing aperture radii as a function of stel- 
lar PSF also counters the problem of variations in measured 
brightness with changing see ing, which is encount ered when 
using small fixed apertures IlLicandro e t al. 2000). In addi- 
tion, photometry of the field stars used to measure the dif- 
ferent ial magnitudes, and of standard stars from the lLandoltl 
( 1992) catalogue, was performed using an aperture of diam- 
eter 10". The smaller aperture was used for differential pho- 
tometry between the comets and field stars, while the larger 
diamete r was used for p hotometric calibration and was that 
used bv lLandoltl Jl992l) . Using the IRAF package PHOT- 
CAL these standard star measurements gave the zero-point, 
extinction coefficient and colour term for each filter for each 
night. These were then used to calculate the magnitudes of 
the field stars in each frame, taking the mean of these val- 
ues gave us a very accurate measurement of the brightness of 
our comparison stars. Adding this value to each of the differ- 
ential comet magnitudes gave us accurate calibrated comet 
magnitudes. There is also a term including the difference in 
colour between the comet and stars involved in the trans- 
formation from differential to calibrated magnitudes, but 
this is negligible when the comparison stars and the comet 
have similar colours. Care was taken to select comparison 
stars with colours in the range expected for cometary nuclei 
to minimise any uncertainties in this colour term, although 
experimentation found this method of obtaining calibrated 
photometry to be robust and independent of which compar- 
ison stars were chosen. 



3 RESULTS FROM SNAP-SHOT IMAGING 

We first obtained 'snap-shot' images of all our potential tar- 
get comets, to ascertain their suitability for time-series pho- 
tometry. For four comets which were not selected as primary 
targets, we were still able to obtain measurements or limits 
on the brightness and activity level of the nuclei from these 
data. These results are summarised in table |2] 

The measured ii-band magnitude tu r is related to the 
effective radius of the nucleus tn (m) by (Russelll ll916t) 

ArtI = 2.238 x 10 22 7?2 A 2 10 o.4(m -™ i? .+,3 Q ) ^ 

where Ar is the geometric albedo and mg = —27.09 is the 
apparent magnitude of the Sun, both in the _R-band. Rh and 
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Table 1. Summary of observations obtained using NTT+EMMI. 



Comet 


UT Date 


n h [A\j] a 


A [AU] 


a [deg.] 


Arf rames X Filter 


Exposure time [s] 


17P/Holmes 


05/03/2005 


4.66 1 


3.93 


9.0 


UxR, 2xV, 2x1 


100 




06/03/2005 


4.66 1 


3.91 


8.9 


25xR, 3xV, 3x1 


100 


43P /Wolf-Harrington 


05/03/2005 


3.30° 


2.61 


13.9 


ixR 


100 


44P/Reinmuth 2 


06/03/2005 


5.17 7 


4.21 


3.0 


2xR 


so 


47P / Ashbrook- Jackson 


05/03/2005 


5.42 7 


4.48 


3.5 


19xiJ, 2xV, 2x1 


85 




06/03/2005 


5.42 1 


4.47 


3.2 


34xR, 2xV, 2x1 


85 


103P/Hartley 2 


06/03/2005 


3.24° 


2.28 


5.0 


3xR, 


45 


104P/Kowal 2 


05/03/2005 


3.06° 


2.31 


13.8 


2xR, 


80 


137P/Shoemaker-Levy 2 


06/03/2005 


6.95 7 


6.17 


5.3 


22xR, 2xV, 2x1 


140 



a Superscripts / and O refer to whether the comet is inbound (pre-perihelion) or outbound (post-perihelion). The variations in 
-Rhi A and a over the course of any one night were smaller than the significance quoted here. 



Table 2. Derived values and limits on radii and activity from snap-shot photometry on four comets. 



Comet 



m R 



m c 771^(1,1,0) rjM [km] Afp [cm] c 



UNDETECTED 

104P 

UNRESOLVED 

44P 

ACTIVE 

43P 
103P 



>23.3° 



>18.. 



<0.56 c 



22.62±0.12 >26.9 15.83±0.12 1.96±0.11 



20.53±0.03 
18.72±0.02 



>15.37 
>14.21 



<2.42 
<4.13 



206±2 
196±1 



a Afp measured through an aperture of radius 5". 
^ 3<x limiting magnitudes. 
c 3a upper limits. 



A are the heliocentric and geocentric distances in AU, and 
a and /3 are the phase angle in degrees and the phase co- 
efficient in magnitudes per degree. Without data acquired 
over a range of phase angles, we must assume a value for 
the phase coefficient. We ta ke the commonly a ssumed value 
of /3 = 0.035 mag. deg" 1 dLamv et alJl2004h . In the ab- 
sence of simultaneous thermal observations, the albedo must 
also be assumed. Comets have been found to be exception- 
ally dark bodies; we take the typical value of Ar = 0.04 
jLamv et al.ll2004l) . We also give the absolute magnitude of 
each nucleus, which allows comparisons without having to 
assume an albedo. This is given by 

m R (l,l,Q)=m R -5log{R h A)-f3a, (2) 

and is the magnitude that would be measured at a hypo- 
thetical point at unit heliocentric and geocentric distances, 
and at a = 0°. 



3.1 Undetected comets 

lLowrv. Fitzsimmons fc Collander-Browr] (120031) have previ- 
ously observed 104P/Kowal 2 when inactive at a heliocen- 
tric distance of 3.9 AU, in snap-shot data acquired with the 
1.0m Jacobus Kapteyn Telescope (JKT) on La Palma in 
June 1999. They obtained a radius measurement of 1.0±0.5 
km, again assuming an albedo of 0.04, which allowed us to 



predict an apparent nuclear magnitude of 21.1 during our 
NTT observations, making this one of our primary targets. 
However, this comet was not immediately detected in our 
initial snap-shot imaging and removed from the target list 
for time-series observations. Upon subsequent processing it 
became clear that the comet was undetectable even in im- 
ages combining both i?-band frames. The comet has a well 
known orbit, having been observed during four apparitions 
since its discovery in 1979 and most recently in February 
2004, and an astrometric fit allowed us to predict the pre- 
cise position of the comet at the time of observation. There 
was no object detected to a ttir = 24.65 ± 0.46 within the 
expected area. Using equation|5|this gives a 3<r limiting mag- 
nitude of mu(l, 1,0) > 18.53, while equation Q allows us to 
put a 3<r upper limit on the radius of the nucleus of 0.56 km, 
assuming Ar = 0.04. 

This limit is in agreement with the previous measure- 
ment, although at the small end of the e xpected size. Comet s 
are known to spontaneously disrupt feoehnhardd I2OO4T) : 
break-up of this nucleus during its perihelion passage in 
May 2004 would also explain our non-detection. However 
Reach et al. observed the comet with the Spitzer infra-red 



1 A Google search found the 104P observations on Yan Fernan- 
dez's useful page of the Spitzer logs for Solar System objects - 
http : //www. physics .ucf . edu/~yf ernandez/sss .html 
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44P/Reinmuth 2 [R] 




Figure 1. Image showing 44P/Reinmuth 2, made up of 2x80 
second exposures. The comet, in the centre of the frame, is very 
faint but appears to be inactive. 




p (arcsec) 

Figure 2. Surface brightness profile of 44P, showing surface 
brightness in magnitudes per square arcsecond against distance 
from the centre of the comet p (arcsec). Beyond 2.5" sky noise 
dominates the profile, as can be seen by the large error bars, how- 
ever the comet (points) and scaled stellar (solid line) profiles are 
indistinguishable in the inner part of the profile, implying that 
the comet was indeed inactive. 



space telescope at almost the same time as our NTT obser- 
vations. The comet was detected at the expected position in 
the Spitzer data (Reach, private communication), so it has 
not broken up and therefore our non-detection gives a true 
upper limit to the size. 

3.2 Unresolved comets 

44P/Reinmuth 2 was expected to have rriR — 22.8 at the 
time of our observations, and was therefore a secondary tar- 
get. Only one frame was taken of this comet on our first 
night of observations before high humidity forced us to close 
the dome, making it impossible to distinguish the comet 
from faint background stars. We acquired two more frames 
on 44P on the second night. 

The comet appears stellar in each of these frames, and 
in a combined frame (Fig. 0. To put limits on any unre- 
solved coma, we measure the surface brightness E at a large 






Figure 3. Images and surface brightness profiles for 43P/Wolf- 
Harrington and 103P/Hartley 2. The comets are visibly active. 
The solid diagonal lines in the top right corner of each profile show 
gradients of -1 and -1.5, from two theoretical models of steady 
state coma; both comets have profiles with gradients ~ — 1. Note 
that the error bars on the surface brightness are smaller than the 
data points. 



distance p (arcsec) from the centre of the comet image. For 
a steady state coma, the surface brightness is inversely pro- 
portional to p, and t he integrated coma magn itude within 
p, m c (p), is given by (IJewitt fc Danielsonl ll984) 



m c (p) = E c (p) - 2.51og(27rp 



(3) 



Figure |5] shows the measured surface brightness profile, 
which appears stellar within the inner 2.5", but is domi- 
nated beyond that by sky noise. We put a limit on any coma 
within this radius at m c (2.5) > 26.9, corresponding to at 
most 2% of the total flux from the comet. We conclude that 
this comet was inactive at the time of observation, which we 
would expect at i?h = 5.2 AU. 

We measure tur = 22.62±0.12, corresponding to 
= 1.96±0.11. This measurement is co nsistent with t he up - 
per limit of r-N < 3.0 km found by lLowrv et all J2003) . 
It is slightly large r than the value of 1.61 km quoted by 
lLamv et al.l i2004T) , but is in agreement if we assume that 
the nucleus is slightly elongated. Such a range in magni- 
tudes would require a minimum axial ratio of a/b = 1.5, 
which is typical of the values determined for JFC nuclei (see 
section |SJ . 



3.3 Active comets 

43P/ Wolf-Harrington and 103P/Hartley 2 were both visibly 
active in snap-shot frames (fig. [3J, and were immediately 
rejected as unsuitable targets for our programme of nucleus 
photometry. We can still place upper limits on the size of 
the nucleus of these comets, as the flux from the nucleus 
cannot be larger than the total flux measured. For 43P, we 
measured a total brightness within an aperture of radius 
5" of m_R = 20.53±0.03, which gives an upper limit on the 
effective radius of tn < 2.42 km. 

The level of ac tivity of a comet can be quantified using 
the quantity Afp jA'Hearn et al.lfl984l) . which is roughly 
proportional to the dust production rate of the comet, as- 
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suming steady state production. It is the product of the 
geometric albedo A, the radius of the photometric aperture 
at the comet p, and the filling factor /, which is the total 
scattering cross-section of the grains within this aperture. 
We calculate this quantity by 



17P/Holmes [R] 



Afp[cm) 



(2AR h ) 2 F ct 



(4) 



where Rh is in AU, but A and p are in cm. F colnc t and 
Fq are the fluxes measured at the Earth from the comet 
are from the Sun respectively. Expressed in this fashion, 
Afp is theoretically independent of the chosen aperture. We 
take an aperture of radius 5", which corresponds to p = 
9400 km for 43P, and gives Afp = 206 ± 2 cm. This is a 
large value, which is expected for a comet which shows high 
activity levels well beyond the canonical 3 AU 'cut-off' for 
H2O sublimation. This comet's perihelion distance jumped 
from q = 2.5 to 1.5 during 1936; it is likely that this major 
change in its orbit is connected to the high activity level 
observed. 

This activity level has frustrated previous efforts to ob- 
tain a definitive measurement of the radius of the nucleus; 
although there have been a number of published snap-shot 
observations of 43P, it has mostly been observed when ac- 
tive, even at distances beyond 4 AU (table |HJ . The upper 
limits obtained whe n the comet was act i ve are all below the 
radii determined by lLowrv et al.l il999L l200St) from images 
in which the comet appeared inactive. There are two possi- 
ble explanations for this discrepancy; firstly, that the range 
in reported magnitudes is due to elongation of the comet, 
in which case the full ra nge of Am = 0.97 mag, from the 
peer-reviewed literature llLowrv et al]ll999t iLicandro et ail 
120001: lLowrv et alJl2003l: lLowrv fc Fitzsimmonsll2005h gives 
a minimum axial ratio of a/b > l0 4Am = 2.4. The alter- 
native explanation is that the apparently stellar observa- 
tions contained significant unresolved coma. This would be 
quite possible, as both were obtained with the 1.0m JKT, 
and the surface brightness profile and equation[3]reveal that 
m c ~ m_R in both cases, suggesting that unresolved coma 
could have contributed over 90% of the flux in each of these 

obs ervations. 

lLowrv fc Fitzsimmonsl <l2005l) argue for the first sce- 
nario, as Lowry's images were taken at larger Rh, the two 
values of rN obtained are consistent with each other and 
a/b > 2.4 falls within the observed range for nuclei (see sec- 
tion EJ. However, our additional limit below 2.5 km pushes 
the weight of probability towards the second explanation. 
Furthermore, if we include Hainaut et al.'s 1996 result of 
my (1 ,1,0) > 16.4 (MPC 27955, qu oted in ILicandro et a" 
2000) and use (V-R) = 0.15±0.08 llLowrv fc Fitzsimmom 



2005) to convert to the _R-band, the first scenario requires 
Am = 1.82, and therefore a/b > 5, which we regard as un- 
realistically large. Even taking the 3<r upper limit on the 
colour of (V ~ R) = 0.39 gives r N < 1.8 km and a/b > 4. We 
therefore conclude that there was significant flux from unre- 
solved coma present in Lowry et al.'s observations, and that 
43P definitely has tn < 2.4 km, and quite possibly rN < 1.6 
km. 

For 103P we found m R = 18.72±0.02, implying r N < 
4.13 km. Within p = 8200 km we measure Afp = 196 ± 
1 cm. Such high activity is unsurprising as the comet was 
outbound at only 3.2 AU, and it has been observed ac- 




Figure 4. Image showing 17P/Holmes, made up of 25x100 sec- 
ond exposures taken on the second night. The comet appears to 
be inactive. The frames are combined in such a way as to remove 
background sky, fixed objects and cosmic rays. 



five at 4.6 AU llLowrv et al]l200St . The radius limit is in 
agree ment with o ther limits of < 5.3 km l|Licandro et ail 
l2000j). < 5.9 km lLowrv fc Fitzsimmonsl ElOlh and < 5^8 
km llLowrv et al.ll2003f) . although it has bee n reported that 
the ra dius is much smaller than these limits; lOroussin et alJ 
(2004) give a radius of 0.71 ± 0.13 km from Infrared Space 
Observatory observations of the comet when it was active 
at R h = 1.2 AU in February 1998. 



4 TIME-SERIES PHOTOMETRY 
4.1 17P/Holmes 

17P/Holmes appears to have been inactive at the time of 
observation. Figure 3] shows a combined image showing only 
the comet, produced by removing a scaled version of the 
stellar background from each frame prior to combining them, 
and is made up of all 25 i?-band exposures taken on the 
second night, giving an effective exposure of 2500 seconds 
= 0.7 hours. A surface brightness profile (fig. |KJ shows no 
sign of activity, and we use equation |!|]to p U t a limit on any 
unresolved coma at m c > 24.6, or < 13.9 ± 28.6 % of the 
total flux. 

We use t he Fourier method described in 
ISnoderass et alJ ll2005t) to search for periodicities in 
the brightness variations of the comet (data given in 
table 0J. We produce a reduced y/ 2 periodogram, where 
minima with reduced x 2 in the range X % l v = 1 ± y^Ju give 
residuals with a Gaussian distribution with la variance, v 
is the number of degrees of freedom of the model, given by 
v = (N ~ 3) where N is the number of data points, for a 
first order Fourier model. 

Figure [§] shows the periodogram for 17P, with x % /v 
against period in d ays. We have a tota l of 39 data points 
in the i?-band, so y/2/u = ^2/(39 - 3) = 0.236. There are 
clearly a number of solutions within the dotted lines marking 
X 2 /v = 1 ±0.236; in fact our null hypothesis (that there is no 
periodic variation, and all scatter is due to noise around the 
mean magnitude) is only rejected at the 2a level. With only 
1.5 nights of data it is impossible to choose between the pos- 
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Table 3. Reported observations of 43P/Wolf-Harrington 



Date 




A 


Q 




m c 


m«(l, 1,0) 6 




App. c 


Ref. 




[AU] 


[AU] 


[deg.] 








[km] 






04/02/1992 


3.04° 


2.69 


18.6 


20.45 ± 0.13 d 




>15.24 


<2.57 


A 


Licandro et al. 2000 


24/08/1995 


4.87 7 


4.26 


10.1 


21.40 ± 0.50 


>21.48 


14.46 


3.67 


S 


Lowry et al. 1999 


June 1996 


3.9 7 










>16.25 e 


<1.61 


A 


Hainaut et al. 1996 


13/06/1999 


4.46° 


3.66 


9.0 


20.81 ± 0.14 


>20.35 


14.43 


3.72 


S 


Lowry et al. 2003 


12/07/2002 


4.43 7 


3.45 


3.8 


21.45 ± 0.06 




>15.40 


<2.38 


A 


Lowry & Fitzsimmons 2005 


05/03/2005 


3.30° 


2.61 


13.9 


20.53 ± 0.03 




>15.37 


<2.42 


A 


this work 



a Superscripts I and O refer to whether the comet is inbound (pre-perihelion) or outbound (post-perihelion). 
^ Recalculated here using rriQ=-27.09, A^=0.04 and /3=0.035, for consistency. 
c Appearance is given as either A — active, or S - stellar. 

d From reported m v = 20.6 ± 0.1, using (V - R) = 0.15 ± 0.08 from Lowry & Fitzsimmons 2005 
e From reported my(l, 1, 0) > 16.4, using (V — R) = 0.15 ± 0.08 from Lowry & Fitzsimmons 2005 




0.5 1 2 5 

p (arcsec) 

Figure 5. Same as fig. |3]for 17P. The stellar profile matches the 
comet, implying that the comet was inactive. 
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Figure 6. Pcriodogram for 17P. Reduced \ 2 is shown against 
period in days. A good fit should fall within x 2 / u '- 
the solid horizontal line is at X 2 / u = lj the dotted line marks the 
lcr level. There are many statistically acceptable minima: those 
marked 'A' are at periods between 30 and 90 minutes and are due 
to the data sampling frequency. The data folded onto the periods 
with each of the four strongest minima are shown in fig. |7| 







Figure 7. 17P photometric data folded onto each of the four 
strongest periods found in the periodogram (fig. These have 
P rot = 7.2 (B), 8.6 (C), 10.3 (D) and 12.8 (E) hours. The 1.5 
nights of data acquired on this comet are insufficient to choose 
one of these periods over the others. 



sible periods. We can reject those marked 'A' in fig.[(j]as they 
are at periods between half and 1.5 hours, and are due to 
the various aliases created by observing the comet in blocks 
with temporal separation of that order. In fig.[7]we show the 
data folded onto the rotation periods corresponding to each 
of the four periods with \ l v ~ 1 (marked B-E), and note 
that they each produce a visually acceptable light-curve. 

We assume that the variation in brightness is due to the 
nucleus being a rotating non-spherical body, and therefore 
expect to see a double peaked light-curve with a rotation 
period P ro t twice the fitted Fourier period Pfltted- The four 
minima are found at Pflttcd = 3.6 (B), 4.3 (C), 5.2 (D) and 
6.4 (E) hours, and correspond to rotation periods of P rot = 
7.2, 8.6, 10.3 and 12.8 hours. 

We plot the data from night 1 as open circles and the 
night 2 data as crosses. While the period remains uncer- 
tain, it can be seen that we observed 17P at a light-curve 
minimum on the first night, and at around maxima on the 
second. Unfortunately there is a gap in the data on the sec- 
ond night, due to the comet passing in front of a star, which 
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appears to have coincided with the minimum between these 
maxima, assuming that one of the above periods is correct. 
We are confident that the calibration of the data from the 
two different nights is accurate, as repeating the method 
with different comparison stars produced no significant dif- 
ference in the calibrated magnitudes; the uncertainty on the 
night-to-night calibration is much smaller than the error 
bars on individual points. Therefore we can be confident 
that we have observed the full amplitude of the light-curve 
at this epoch. Here we follow the common practice of de- 
scribing the nucleus as the simplest non-spherical shape: a 
tri-axial ellipsoid with semi-axes a > b — c. The variations in 
brightness that produce the light-curve are then due to the 
changing cross-section of the rotating elongated nucleus. We 
can only measure a lower limit on a/b, as we do not know the 
orientation of the rotation axis. From Am = 0.3 we obtain 
a lower limit on the elongation of the nucleus of a/b > 1.3. 

We measure the mean apparent magnitude of 17P to be 
ttir = 22.86 ± 0.02. Using equation^ this gives an effective 
radius of tn = 1.62 ± 0.01 km for the equivalent spherical 
body. Taking a/b = 1.3, this corresponds to dimensions of 
the nucleus of a x b — 1.9 x 1.4 km. This is in excellent 
agreement with the tn = 1.71 km measurement fr om a snap- 
shot observation acquired bv iLamv et alJ <l2000Tl using the 
HST. 

Measuring the elongation from the range of the light- 
curve requires that the variations are in fact due to changes 
in the projected surface area of the comet, and not due to 
large scale changes in albedo. To test this assumption, we 
measured the colour indices (V — R) and (R—I) at least once 
per night, to check for any colour variations which would 
indicate varying surface properties. We acquired 'colour 
blocks' by cycling through the three filters VR and / in the 
following order: RVRIR, repeating the last four frames for 
faint objects. This allowed direct interpolation to find the R- 
band magnitude at the time of the V or /-band observations, 
by simply averaging the i?-band frames bracketing them. For 
17P we acquired one colour block, containing 2 x 100s V and 
/-band frames in the sequence RVR1 RVRIR, on the first 
night and two blocks on the second. The colours measured 
were consistent within the uncertainties on each, supporting 
our assumption that variations are due to the shape of the 
nucleus, and give average colours of (V — R) — 0.41 ± 0.07 
and (/?-/) = 0.44 ± 0.08. 

4.2 47P/Ashbrook- Jackson 

47P/Ashbrook- Jackson was visible throughout the entire 
time we were able to observe on both nights, and conse- 
quently was well monitored, with a total of 53 _R-band ex- 
posures of 85 seconds each (table |SJ . The comet appeared 
bright and stellar in each; a combination of all of the frames 
taken on the second night is shown in fig. HI A surface bright- 
ness profile (fig.|5Jl shows that the comet has negligible coma, 
with a scaled stellar profile matching the comet down to sur- 
face brightness £ > 30 mag arcsec -2 . Using equation [3] we 
put limits on any unresolved coma at m c > 27.4, imply- 
ing that activity is entirely negligible and can contribute no 
more than 0.5 ± 2.1% of the total flux, as was expected as 
this comet was close to aphelion at the time of observation. 

There are four statistically acceptable minima in the pe- 
riodogram for 47P, which are marked A-D in fig. 1101 Light- 



47P/Ashbrook-Jackson [R] 




Figure 8. Image showing 47P/Ashbrook- Jackson, made up of 
34x85 second exposures taken on the second night. The comet 
appears to be inactive. 
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Figure 9. Same as fig.|2]for 47P. The stellar profile matches the 
comet, implying that the comet was inactive. 
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Figure 10. Periodogram for 47P. There are four statistically ac- 
ceptable minima: the data folded onto the periods with each of 
these minima are shown in fig. 1111 
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Figure 11. 47P photometric data folded onto each of the four 
strongest periods found in the pcriodogram 

(fig-Eni- These have 

P rot = 11.2 (A), 15.5 (B), 21.6 (C) and 44.0 (D) hours. The 1.5 
nights of data acquired on this comet are insufficient to choose 
one of these periods over the others. 



curves showing the data folded onto each of these are shown 
in fig. 1111 Period 'C (P ro t = 21.6 hrs) is marginally the 
strongest minima, but the data do not allow us to choose 
between these periods. In this figure it can clearly be seen 
that over the course of the second night a definite maximum 
was observed, and the approach to a maximum appears to 
have been observed on the first night, which puts strong con- 
straints on the rotation period, although we cannot choose 
between the four possible periods shown based on these data. 
Note that the data do allow us to measure the period pre- 
cisely, with small error bars on the determined periods of 
11.2 ± 0.3, 15.5 ± 0.5, 21.6 ± 1.0 and 44.0^2:9 hours > but do 
not give a unique solution. It is clear that the periods are 
aliases of each other at ~ 1, 3/2, 2 and 4x11 hours, yet it is 
not possible to select one as the correct period. The rotation 
period cannot be longer than 44 hours (period 'D'), as this 
period requires that we saw alternate maxima of a double 
peaked light-curve on the 2 nights; the remaining solutions 
imply that there were other maxima unobserved between 
these events. There is an unpublished partial light-curve on 
47P acquired by Lamy et al. using the HST, which puts a 
lower limit on the rotation period at 44.5 hours. If this is 
correct, it would suggest that period 'D' is the real one, and 
that the shape of the body is irregular: The 'best-fit' sinu- 
soid drawn through the data to guide the eye in fig. IllT D) 
does not represent the data well. The range in magnitudes 
from the data is Am = 0.45, giving a/b > 1.5, which is 
typical of JFC nuclei. 

The mean magnitude was measured to be tur = 21.68± 
0.01, which gives a radius of 3.38 ± 0.01 km, assuming an 
albedo of 4%. Taken with the elongation implied by the full 
range in the data, a/b > 1.5, we get dimensions of the nu- 
cleus of a x b — 4.2 x 2.8 km, for what is most likely an 
overly simple model of the nucleus. T his is in agreemen t 
with the radius measurements of 2.8 km llLamv et alj|2004l) . 



3.1 km (|lJ£androe^alJ 
found by lLowr^Tt^uT 

(V-R) 



200; 

0.45 ± 0.03 and (R 



and the limit of tn < 6.1 km 
. The colours measured were 
/) = 0.38 ±0.03. 



Figure 12. Image showing 137P/Shoemaker-Levy 2, created 
from 22x140 second exposures taken on the second night. 
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Figure 13. Same as fig. [5] for 137P. The stellar profile matches 
the comet, implying that the comet was inactive. 



4.3 137P/Shoemaker-Levy 2 

137P/Shoemaker-Levy 2 was a secondary target, as its ex- 
tremely large distance (7.0 AU, inbound) gave it a predicted 
brightness of ttir — 23.0, making it a challenging target even 
with the 3.6m NTT. However, we attempted a snap-shot on 
the second night of observations at a point when all our pri- 
mary targets were either below the horizon or too close to 
field stars. The comet was detected at a similar brightness to 
17P, so was added to our light-curve target list. An image 
produced by shifting and adding all i?-band data and the 
corresponding profile (figs. 1121 & 1131 show that the comet 
was inactive, as expected as it was close to aphelion. The 
limit on any coma contribution within a 5" aperture is for- 
mally found to be m c > 22.9, which corresponds to up to 
82% of the flux, however at p = 5" sky noise entirely dom- 
inates; although most of the background sky is removed by 
the process that gives fig. 1121 the comet profile can not be 
measured accurately beyond £ = 32 mag arcsec -2 . The pro- 
file is clearly stellar in the inner part; measuring the surface 
brightness at p — 2" gives m c > 25.0, corresponding to 
< 14± 10% of the flux. 

We took 22 i?-band images of 137P, however the loss of 
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Figure 14. Photometric data for 137P. There is no statistically 
significant variation within the data; in the one night of data we 
find no evidence for rotation. 



the 3rd night to bad weather meant that this data covered 
only 3 hours (table There is no statistically significant 
variation within this data, therefore we cannot constrain the 
rotation period. We show the data in fig. 1141 This time- 
series allows us to obtain an accurate measurement of the 
brightness of this inactive nucleus, and therefore a size, again 
assuming a typical 4% albedo. We measure vtir — 22.86 ± 
0.03, and therefore tn = 3.58 ± 0.05 k m, which is similar t o 
the limit of < 3.4 km measured bv lLowrv et alJ {2003). 
We also measured colour indices for this nucleus of (V — 
R) = 0.71 ± 0.18 and (R - I) = 0.54 ± 0.15, which places 
137P at the red end of the distribution of measured nuclei 
colours, admittedly with relatively large error bars due to 
its faintness. 



5 DISCUSSION OF COMET PROPERTIES 
AND COMPARISON WITH KBOS 

5.1 Bulk densities 

We can use the rotation period and elongation found for 
17P and 47P to put limits on the bulk density of the nu- 
cleus Dn (g cm -3 ). By assuming that nuclei have neg- 
ligible tensile strength [as evidenced by the break up of 
D /Shoe maker- Levy 9 under th e gravitational influence of 
Jupiter dAsphaug fc Bend ^e) and results from the Dee p 
Impact mission to comet 9P/Temple 1 llA'Hearn et all 2005 ) ; 
see also IWeissman. Asphaug fc Lowrvl J2004f) l a lower limit 
on the density can be found by balancing self gravity 
with centrifugal force. We use the approximation given by 
IPravec fc Harrisl (l2000l) ; 



D N > 



10.9 a 
P 2 b 

-'rot u 



(5) 



where P ro t is in hours. This equation assumes that the grav- 
itational acceleration at the end of the long axis a of the 
ellipsoid is the same as that for a sphere, but reduced by 
a factor of b/a. It provides an acceptable approximation to 
the more complex s olution found by integrating over the vol- 
ume of an ellipsoid (iRichardso n. Elankumaran fc Sanderson! 
2005) 
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Figure 15. Rotation period against projected axial ratio for all 
JFC nuclei. The error bars on 17P and 47P cover the full range 
of possible periods discussed in the text; the points shown are 
for the nominally 'best' periods. The solid lines show where a 
strengthless body would lie on these axes for a variety of densities. 
There appea rs to be a cut-off in densit ies at ~ 0.6 g cm -3 , as 
sugg ested bv [Lo wry & ^yeissmar] j 2003). Data on other comets 
from lSnodgrass et al.ll2005t and references therein. 
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where G is the gravitational constant and e = y/l — (^) 2 . 
T he derivation of thi s equa tion i s given i n two different ways 
bv IRichardson et all (120051) and lHarrisI (1200% it is also en- 
tirely equivalent to the expression given by iLuu fc JewittJ 
lll992F) [their equation 9]. 

Using equation |5] we can calculate the minimum bulk 
density of each nucleus using the light-curve derived param- 
eters a/b and P ro t- The calculated densities, along with all 
other results from the time-series data, are presented in ta- 
ble Q The full range of possible periods give densities Dn > 
0.09 g cm" 3 for 17P, and D N > 0.01 g cm" 3 for 47P. Even 
if spinning at the shortest of the possible periods found for 
these comets, they do not require densities larger than 0.28 
and 0.13 g cm" 3 . All of these lower density limits are con- 
sistent with the range observed for other nuclei; in this re- 
spect the nuclei of 17P and 47P are fairly typical, for any 
of the possible values of their rotation periods. The rela- 
tively long periods do not put strong constraints on their 
individual densities, however considering all measurements 
obtained we can infer valuable information concerning the 
densities of the JFC population. Under the assumption that 
all nuclei have a similar Dn, then the values obtained by 
the light-curve technique will generally be below this value, 
but none will be observed above it. For a large enough sam- 
ple of a population with randomly inclined rotation axes 
and a spin rate distribution, any observed upper limit to 
the minimum densities derived will give the true density of 
the nuclei. Figure 1151 shows the logarithm of the rotation 
period against projected axial ratio for all JFC nuclei. The 
plot shows that there appears to be a cut-off in densities at 
around 0.6 g cm" 3 . This agrees with the results from the 
Deep Impact mission, which found c omet 9P/Temple 1 to 
have a low density of 0.4±0.3 g cm" 3 (IRichardson fc Meloshl 
2006). In asteroids, there is a clear cut-off for bodies with 
r > 75 m at ~ 3 g cm " 3 (IPravec et al]l2002l) . and it has 
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Table 7. Derived physical parameters and colours from time series photometry on three comets. 



Comet 


m R 


m c 


m fl (l,l,0) 


r N 
[km] 


Prot a 

[hr] 


a/b b 


[g cm 3 ] 


(V-R) 


(R-I) 


17P 
47P 
137P 


22.864±0.020 
21.680±0.007 
22.862±0.032 


>24.6 
>27.4 
>22.9 


16.241±0.020 
14.638±0.007 
14.513±0.032 


1.62±0.01 
3.38±0.01 
3.58±0.05 


7.2—12.8 
11.2—44.0 

n/a 


1.3±0.1 
1.5±0.1 

n/a 


0.09±0.02 
0.01±0.01 

n/a 


0.41±0.07 
0.45±0.03 
0.71±0.18 


0.44±0.08 
0.38±0.03 
0.54±0.15 



The rotation periods of both 17P and 47P are found to be one of four discrete values within the given ranges (see text). 
Lower limits as the orientation of the rotation axis is unknown. 



been suggeste d l|LowrW^ye issmanll2003l : IWeissman et alJ 
l2004l:ISnodgrass et al.ll2005h that the 0.6 g cm" 3 cut-off for 
comets is equivalent to this, the lower value reflecting on the 
fact that comets are under-dense when compared to aster- 
oids, due to their larger volatile content and/or more porous 
structure than asteroids. 

The relative proportions of ice and rock within a body 
and the porosity cannot be found independently, however 
we can express the fraction of a body, by mass, made up of 
ice in terms of Dn and the porosity tp, which is the fraction 
of void space within the total volume: 

M icc _ 1 - (1 - tp)D Iock /D N 



Ah 



1 — -Drock/Mc 



(7) 



where D lcc and -Drock refer to the densities of the volatile ices 
(~ 1 g cm -3 ) and silicates (~ 3 g cm -3 ) respectively. Using 
this equation for a range of assumed values of tp allows us to 
calculate the corresponding proportions of ice and rock that 
would be required to give a density of 0.6 g cm -3 . Inserting 
this value for Dn into equation |7] gives the mass fraction of 
ice of Mi C e/MN = 2 — 2.5tp, and thus a minimum value for 
the porosity of tp = 40% for an entirely ice nucleus. This 
is considerably higher than the values found through lab- 
based studies of meteorites, which have an average micro- 
poros ity around 10%, with the m ost porous having tp w 
30% llBritt fc Consolmagnoi r2003b suggesting considerable 
macro-porosity or voids within nuclei. A value of ip = 50% 
requires that the nucleus is 75% ice to give Dn = 0.6 g cm -3 , 
while a nucleus with this density and a porosity of ip > 60% 
must be silicate dominated (less than half of its mass as ice) . 
In the limiting case of a comet almost entirely depleted in 
volatiles (Mi ce — + 0), a ~ 100% silicate nucleus must have 
ip ~ 80% to give the low observed density. 

There is one object which has a minimum density 
greater than 0.6 g cm~ 3 in fig.HSI 133P/Elst-Pizarro. The 
presence of 133P does not present a strong case against 
the 0.6 g cm -3 cut-off, as this object is thought to be an 
asteroid-comet transition object; it was first found and num- 
bered as asteroid 7968, and has an asteroid-like orbit, but 
has subsequently been found to have a dust trail along its or- 
bit, and shows sp oradic activity llHsieh. Jewitt fc Fernandez! 
l2004l;lTothll200fT) . This and similar objects have been shown 
to be from a dynamically and p ossibly chemically distinct 
population llHsieh fc Jewittl2006t) . Their stable orbits imply 
that these Main Belt comets have remained at or near their 
formation site. Due to the observed compositional gradient 
in the Solar System, this would imply that they possess a sig- 
nificantly higher silicate content than the JFCs that formed 
in the Kuiper Belt, and hence a higher density. Therefore 
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Figure 16. Same as fig. 1151 but includin g data on KBO s (fille d 
circles) and Centaurs (tri angl es') from [ Lacerda fc Luul feOOrJ) . 
iTrilling fc Bernstein! <2006l) and lOrtiz et alJfcOOd) . The distribu- 



tion of JFCs is very similar to that of KBOs/Ccntaurs, although 
the KBOs tend to have lower axial ratios, most likely due to the 
larger size of known KBOs. The KBO data also falls above the 
0.6 g cm -3 density line, and taken together with the JFC data 
support the hypothesis of a cut-off at this density. The labelled 
objects are referred to in the text. 



we disregard 133P when discussing the global properties of 
he JFC population alone. 

We compare the position of JFC nuclei in fig. 1151 with 
the results from similar studies on their supposed parent 
population, the KBOs. We select KBOs with Hv > 5.1 
(corresponding to bodies with radius < 200 km assuming 
a 10% albedo typical of KBOs) from the recently publ ished 
compilations of all KBO da ta bv lLacerda fc Luul (I2006D and 
ITrilling fc Bernstein! l|2006h [see these papers for references 
to results for i ndivid ual objects], and also the new results 
of lOrtiz et all J20061) . Figure ITBI shows the position of the 
16 KBOs (filled circles) and 14 Centaurs & Scattered disk 
objects (triangles) with measured light-curves from these 
papers together with the JFC nuclei from fig. 1151 The ro- 
tation periods plotted assume that the changes in bright- 
ness are due to the observed cross-section of a rotating elon- 
gated body (i.e. double peaked light-curves). Furthermore, 
to allow comparison between the nuclei and the KBOS, we 
must assume that the KBOs are also effectively strengthless. 
Whether this assumption can be applied to such large bod- 
ies is a n interesting question, however the work of lDavidssonl 
( 1999) shows that for similar rotation periods, larger bodies 
will fall into the 'damaged region' which have been frac- 
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tured due to sheer forces and are therefore only gravitation- 
ally bound. For the very largest bodies there is likely to be 
considerable non-homogeniety due to high central pressures. 
For this reason we restrict the KBOs to the smaller ones for 
which light-curves are available. 

It is clear from fig. 1161 that the comets and the outer 
solar system bodies populate a similar region of axial-ratio- 
period space, which implies that they have similar bulk den- 
sities if it is assumed that neither have significant strength, 
supporting the theory that JFCs come from the Kuiper Belt. 
We note that the KBOs also have minimum bulk densities 
below the postulated 0.6 g cm -3 cut-off; the fastest rotating 
KBO known is 2001 CZ31, which has a rotation period of 
frot = 4.7 hours jLacerda fc Luul2006l) . and therefore a den- 
sity of Dn > 0.60 g cm~ 3 . This agrees with the theory that 
KBOs and comets are rubble piles and under-dense com- 
pared with asteroids, again suggesting a higher proportion 
of volatile ices and higher porosity. 



5.2 Rotational properties 

We now compare the distributions of comet and KBO axial- 
ratios and spin rates separately, to look for similarities 
in shape and in collisional evolution. KBOs also tend to 
have low elongations, generally with minimum axial ratios 
a/b < 2. The exceptions are 2003 BF91 and 2001 QG298, 
which are on the right of the plot with a/b > 2.73 and 
2.86 respectively. 2003 BF91 is one of the smallest KBOs 
known, with an absolute visual magnitude of Hy = 11.7, 
corresponding to a radius of 9 km assuming a 10% albedo 
(Trilli ng fc Bernsteinll2006l) . or r = 15 km assuming a comet 
like 4% albedo. This small size and large axial ratio sug- 
gests that there could be more elongated objects among the 
smaller KBOs, which have not been detected due to the ex- 
treme faintness of these objects. 2001 QG298 also shows a 
large photometric variation, which would correspond to a 
very large elongation for a single r ~ 120 km body, how- 
ever this object is thoug ht to be a contact binary due to the 
shape of its light-curve JSheppard fc Jewittlliool) . 

Figure 1171 shows histograms of spin rates for JFCs, 
KBOs and Centaurs. We investigate the rotation rate dis- 
tribution of JFCs to determine whether they are a colli- 
sionally relaxed population, or whether torques due to out- 
gassing substantially alter this distribution. For a collision- 
ally relaxed population, these distributions should show a 
Maxwellian shape, as seen in large (r > 20 km) asteroids 
iPravec et al 1 l2002l - t heir fig. 3 is re c reated in fig |17f D'l for 
comparison). Following iPravec et al.l (|2002), we plot the ro- 
tation frequency / = 1 /P ro t normalised using the geometric 
mean (/) of each sample, which accounts for different sized 
bodies. 

The histograms are not obviously Maxwellian in nature, 
in fact the comet distribution does not significantly differ 
from a flat distribution; a Kolmogorov-Smirnov test compar- 
ing the normalised frequency distribution in fig. I17f A) with 
a flat distribution with the same mean gives a D statistic of 
0.12, implying only a 6% probability that the distribution is 
not flat. Naively, we expect that torques due to out-gassing 
from nuclei will act to flatten the distribution from an ini- 
tial Maxwellian, as torques will speed up or slow down the 
rotation, although detailed numerical modelling is required 
to give the theoretical distribution for a cometary popula- 




Figure 17. Histogram of normalised spin rates (//{/), where 
/ = l/P rot ) rates for (A) JFCs, (B) KBOs, (C) Cen taurs and (D) 
astero ids. The asteroidal data is reproduced from IPravec et alj 
12002), and shows the Maxwellian distribution in spin rates ex- 
pected for a collisionally evolved system. 
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Figure 18. (R — I) against (V — R) for all JFCs with known 
colours. Filled circles — this work; open circles - comets with 
colours previously determi ned using the same multi-filter pho- 
tometry method used here (Snodgrass et al. 2005, and references 
therein). The position of the Sun on these axes is marked using 
the symbol 0. 



tion. Given the small sample size, further light-curves are 
required for cometary nuclei. 



5.3 Surface colours 

As described in section [I] we take multi-filter images of the 
comets at different epochs to test for variation in colours 
which would imply albedo variations over the surface. We 
find that all of our comets have consistent colours, implying 
that it is shape, and not large scale albedo features, that is 
responsible for the observed brightness variations. The mea- 
sured (V — R) and (R — I) colours are given in table [7| and 
plotted in fig. 1181 These data fit amongst the distribution 
of known nuclei, with the exception of 137P, which is con- 
siderably more red in (V — R) than other comets, although 
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• • •• • 
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(V-R) 



0.8 



Figure 19. (fl - /) against (V - R) for all JFCs, KBOs and 
Centaurs with known colours. Open circles - JFCs; filled circles - 
KBOs; triangles — Centaurs. The position of the Sun on these axes 
is marked using the symbol 0. The dashed line shows our best fit 
to the comet data, the dotted line i s the fi t to the KBOs. KBO an d 
Centaur data from ljewitt fc Luul <200ll) : Ipeixinho et al] <20u4) . 



with relatively large error bars which bring it into touch 
with the others. All three of the comets with new colours 
presen ted here fit along the trend found bv lSnoderass et ail 
(2005), which describes an increase in (R — I) with in- 
creasing (V — R) indicative of a steadily increasing spectra 
through these bands, as seen in primitive asteroids and the 
few c omet ary nu clei which have been observed spectroscopi- 
cally (L uu 1993). The addition of the new nuclei gives a best 
fit straight line to the observed trend: 



(R ~ /)„uc = m(V - R)nu 



m = 0.94 ±0.16 , . 

' ^ c = 0.01 ±0.01 IM 



which is identical to the fit given bv lSnoderass et ail |2005), 
with a very slightly reduced uncertainty on the gradient m. 
In both of these fits, asteroid-comet transition object 133P 
and comet 45P/Honda-Mrkos-Pajdusakova are rejected as 
outliers, as they have (R — I) < (V — R), indicating a down- 
turn in the spectra typical of stony asteroids. 

Figure 1191 is a colour-colour plot showing the data 
on KBOs (filled circles) and that of JFCs (open circles), 
with the same trend line over-plotted (KBO data from 
IJewitt fc Luul200 ltlPeixinho et al.l2004ll . It is clear that this 
line provides a reasonable fit to the data, although there is 
considerable scatter in the KBO results, which is confirmed 
by a fit to the KBO data which is nearly identical to that 
for JFCs: 



(R - J)kbo = m(V - R)kbo + < 



m - 
c = 



0.94 ± 0.06 
-0.03 ±0.01 



(9) 



This is plotted as a dotted line in fig. 1191 However, while 
both populations follow the same trend in colours, it is clear 
that the KBOs are generally redder than the comet nuclei. 
The mean colour of JFC nuclei is (V - R) nuc = 0.46 ± 0.12 
(N = 23), compared with (V - R) KBO = 0.59 ± 0.08 (N = 
62) for the KBO population. 

It has been suggested that 'space weathering' produces 
the very red surfaces seen in KBOs and that cometary ac- 
tivity resurfa ces nuclei to give them more neutral colours 
Jjewitdl2002l . Our analysis supports this hypothesis; the 



similar (V — R)-(R — /) relations suggesting similar spec- 
tra and therefore a common source, while the JFCs tend 
towards the blue side of the distribution and lack the 'ultra- 
red matter' of KBOs. Under this paradigm, the very red 
colour we report for 137P may indicate that it is a rela- 
tively young comet which has undergone little resurfacing, 
although we remain cautious of inferring too much from the 
colour in light of its large uncertainty. We note that the 
Centaurs (triangles in fig. 1191 data also from ljewitt fc Luul 
1200 ltlPeixinho et al.ll2004l) have a similar spread in (V - R) 
to the KBOs, but do not appear to show the same trend 
as the JFCs and KBOs, as they exhibit a very large vari- 
ation in (R — I). This is somewhat surprising, as Centaurs 
are thought to be the intermediate dynamical population 
through which KBOs pass to become JFCs, but there are 
only nine Centaurs in fig. 1191 and further data may clarify 
the situation. 

Given the evidence for the loss of extremely red matter 
due to sublimation and resurfacing, it is important to search 
for observational signatures of continued s urface ev o lution 
within the JFC population. For example, iHughesI <|2003T) 
claimed that the size distribution index was shallower for 
comets with a perihelion distance q > 2.7 AU. If redden- 
ing of surface material occurs on a time-scale shorter than 
that for significant orbital change of a JFC, then one might 
see a correlation with orbital parameters. However, plots of 
(V ~ R) against both semi-major axis a and the fraction of 
orbital period spent at Rh < 3 AU showed no correlation. 
Therefore the nucleus colour does not depend on the orbital 
parameters of the JFC. 

A remaining possibility is that the timescale for signif- 
icant reddening of the surface is shorter than the orbital 
period. 28P/Neujmi n 1 has a large amoun t of extant pho- 
tome try (collated bv lDelahodde et al]l200l] and lLamv et al] 
|2004|) and is well known as a comet that is essentially dor- 
mant throughout much of its orbit, only showing a small 
amount of activity near perihelion. Taking these colours and 
plotting them against position in the orbit shows no corre- 
lation, it appears that the surface colour of 28P is constant 
throughout its orbital cycle. Therefore we conclude that sig- 
nificant reddening does not occur on time-scales of less than 
20 years and the colours measured for individual comets are 
not dependant on their orbital position at the time of ob- 
servation. 



6 SUMMARY 

From time-series VRI photometry of comets 17P, 47P and 
137P, acquired during March 2005 using the 3.6m NTT at La 
Silla, Chile, we place constraints on the sizes, shapes, bulk 
densities and surface colours of their nuclei. Additionally we 
present results from snap-shot observations of comets 43P, 
44P, 103P and 104P. We find: 

(i) Comets 43P and 103P were observed to be active. We 
place upper limits on the radii of their nuclei at tn < 2.4 km 
and < 4.1 km respectively. 104P was not detected, allowing 
us to put a 3cr upper limit of its radius of rN < 0.56 km, 
assuming an albedo of 4%. 

(ii) The remaining detected comets were stellar in appear- 
ance. Surface brightness profiles confirmed that all were ef- 
fectively inactive. By assuming an albedo of 4% and a phase 
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coefficient of /3 = 0.035 mag deg , we convert our mean 
measured magnitudes into radii of 1.96 ± 0.11 km for 44P, 

I. 62±0.01 km for 17P, 3.38±0.01 km for 47P and 3.58±0.05 
km for 137P. 

(iii) From our data, we cannot uniquely determine the ro- 
tation periods for 17P and 47P, but can constrain the values 
to a number of solutions. For 17P, the period lies in the re- 
gion 7.2 < P ro t < 12.8 hours, for 47P the period is within 

II. 2 < P r ot < 44.0 hours. We measure minimum axial ratios 
of a/b > 1.3 and 1.5 for these two comets, implying mini- 
mum densities of between 0.09 and 0.28 g cm -3 for 17P, and 
0.01 to 0.13 g cm" 3 for 47P. 

(iv) KBOs have similar rotation periods and elongations 
to the JFC sample, implying that they have similar densities. 
Both obey a cut-off in bulk density near 0.6 g cm . The 
limited amount of data at present is also consistent with 
them having the same spin-rate distribution. These physical 
constraints support the theory that JFCs come from the 
Kuiper Belt. 

(v) We measure colour indices for four nuclei, three of 
which are typical of other JFCs; (V - R) = 0.41 ± 0.07 
and {R-I) = 0.44 ± 0.08 for 17P, (V - R) = 0.45 ± 0.03 
and {R-I) = 0.38 ± 0.03 for 47P. For 137P, we measure 
(V - R) = 0.71 ± 0.18 and {R-I) = 0.54 ± 0.15, which 
give this comet the reddest nuclear colours measured, albeit 
with large uncertainties. 

(vi) JFC colours follow a trend of increasing {R — I) with 
increasing {V — R) that suggests a continually rising spec- 
trum through these bands. This trend also fits KBOs, imply- 
ing that their spectra have a similar shape. Hence spectrally 
the JFC population is entirely consistent with being derived 
from the KBO population but lacking in the extremely red 
objects. The tendency of nuclei towards the blue end of this 
trend shows that they are depleted in the 'ultra-red matter' 
that is built up on bodies in the outer solar system by 'space 
weathering'. 

(vii) The optical colours of JFCs show no correlation with 
orbital parameters. We also note that the extensive photom- 
etry of 28P/Neujmin 1 in the literature shows no significant 
variation in colour as a function of time since it was previ- 
ously active. Therefore the diversity of JFC nuclei colours 
is not due to either the orbit or to the time at which the 
comet was observed. 
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Table 4. Apparent _R-band magnitudes for comet 17P. 



Table 5. Apparent i?-band magnitudes for comet 47P. 
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Table 6. Apparent _R-band magnitudes for comet 137P. 
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